Walking near a Conformal Fixed Point: the 2-d 0(3) Model at 9 ~ it as a Test Case 
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Slowly walking technicolor models provide a mechanism for electroweak symmetry breaking whose 
nonperturbative lattice investigation is rather challenging. Here we demonstrate walking near a con- 
formal fixed point considering the 2-d lattice 0(3) model at vacuum angle 9 ta ir. The essential 
features of walking technicolor models are shared by this toy model and can be accurately inves- 
tigated by numerical simulations. We show results for the running coupling and the beta-function 
and we perform a finite size scaling analysis of the massgap close to the conformal point. 



The scalar Higgs field that drives spontaneous elec- 
troweak symmetry breaking in the Standard Model is 
considered unnatural as a fundamental degree of free- 
dom because it suffers from the gauge hierarchy prob- 
lem. Technicolor provides a promising mechanism that 
stabilizes the electroweak scale against the Planck or 
GUT scale by introducing a new asymptotically free 
strongly coupled gauge theory The chiral conden- 
sate of techniquarks then induces electroweak symme- 
try breaking and replaces the fundamental Higgs field 
in a natural way, i.e. without fine-tuning. The origi- 
nal technicolor models [2-4] suffer from flavor-changing 
neutral currents. In addition, electroweak precision tests 
are strongly affected by the running of the technicolor 
gauge coupling 5]. In order to avoid these problems, 
it has been proposed to consider slowly walking tech- 
nicolor theories In these models the running of 
the gauge coupling slows down due to the proximity to 
a conformal fixed point Besides non-Abelian tech- 
nicolor gauge fields, these models usually contain many 
flavors of techniquarks in the fundamental representation 
or tcchnifcrmions in higher-dimensional representations 
of the technicolor gauge group. In order to naturally 
stabilize the electroweak scale against the Planck scale, 
these models are still asymptotically free. They are just 
outside the conformal window of theories that have an 
infrared conformal and an ultraviolet asymptotically free 
fixed point. The lower edge of the conformal window 
corresponds to the merging of two fixed points and thus 
to a double zero of the (3- function 1(| 11|- Determining 
whether a theory is inside or near the conformal win- 
dow requires nonperturbative investigations, which can 
be performed from first principles using lattice gauge 
theory. Lattice theories with different gauge groups and 
with different fermionic matter content have been stud- 



gating whether a theory is inside the conformal window 
is a difficult task. Determining whether a theory is slowly 
walking is even more challenging. At present, there is no 
consensus on which theories are slowly walking. 

In this paper, for the first time we unambiguously 
demonstrate walking near a conformal fixed point in an 
asymptotically free field theory — the 2-d 0(3) model, 
which shares many features with 4-d non-Abelian gauge 
theories. In particular, besides being asymptotically free, 
it has a dynamically generated massgap, instantons, and 
thus a non-trivial vacuum angle 9. At 8 = 0, the Eu- 
clidean action of the model is given by 



d 2 x d,,e- d„e. 



(1) 



ied quite intensively |12H18l |. Due to finite-size and finite 
lattice spacing effects as well as due to difficulties with 
simulating nearly massless dynamical fermions, investi- 



Here e(x) is a 3-component unit-vector field and g is the 
coupling constant. The topological charge, 

Q[e\ = ^Jd 2 x e„„e- (9,ex d v e) e U 2 [S 2 } = Z, (2) 

contributes an additional term i9Q[e\ to the action. The 
2-d O (3) model is integrable at 9 — 1^-21 as well as at 
9 = t: [22J, but not at intermediate values of 9. The mass- 
gap at 9 = has been determined analytically [23| and 
is given by M = fA^rg, where e is the base of the nat- 
ural logarithm and A^j- is the scale that is dynamically 
generated by dimensional transmutation in the modified 
minimal subtraction scheme. As one varies 9, the mass- 
gap is reduced until it finally vanishes at 9 = n. The 
conjectured exact S-matrix at 9 = it has recently been 
confirmed by lattice simulations with per mille level accu- 
racy 24|. This study has also demonstrated beyond any 
reasonable doubt that 9 is a relevant parameter that does 
not get renormalized nonperturbatively. Consequently, 
each value of 9 characterizes a different physical theory, 
which was further supported by [25| . 

As one infers from the exact S-matrix [2^, at energies 
far below A-^s the 2-d 0(3) model at 9 — tt reduces to 
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FIG. 1. Triangular lattice of size L x ft. The triangles {xyz) 

carry the topological term i6q^ xyz ) . The correlation function 

C(ti,t2\8) is determined between the time-slices at ti andti. 



the k = 1 WZNW model [26|-[28|. a conformal field theory 
with the Euclidean action 

S[t/] = i f d 2 xTr[d^d fl U}-27rikS WZNW [U}. (3) 
2.9' J 

Here U (x) £ SU(2) — S 3 and g' is a coupling constant. 
The Wess-Zumino-Novikov-Witten term is given by 



Swznw[U 

1 
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d 2 x dx 3 e^pTrp^daUU^d^UU^dpU]. (4) 



Here H 3 is a 3-d hemisphere whose boundary dH 3 = S 2 
is the compactified 2-d space-time. Since ^[S 13 ] = {0}, 
there are no topological obstructions against extending 
the 2-d field U(x) on x £ S 2 to a 3-d field U(x,x 3 ) on 
(x,x 3 ) £ if 3 . The WZNW-terms corresponding to two 
different extensions and differ by an integer, 



Swznw[U^} — Swznw 
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Tr[U^du,UU^d u UU^d p U] e n 3 [S 3 ] = Z. 



(5) 



Here two hemispheres have been combined to form a 
sphere S 3 with U corresponding to on one and to 
on the other hemisphere. Since Swznw[U] is thus 
well-defined only up to an integer, in order to obtain an 
unambiguous contribution exp(2irikSwzNw\U]) to the 
functional integral, the level k must be quantized in inte- 
ger units. Interestingly, the WZNW model has a global 
SU(2) L x SU(2) R symmetry, U(x)' = LU(x)R\ which 
extends the 0(3) symmetry to 0(4). It should be noted 
that only the low-energy WZNW sector, but not the en- 
tire 0(3) model at 9 = tt has the enlarged 0(4) sym- 
metry. The scale A^j^, which results from anomalous 
scale breaking and dimensional transmutation, still in- 
duces explicit symmetry breaking down to 0(3). Due to 
the enlarged 0(4) symmetry of the low-energy sector, an 
0(3) singlet becomes degenerate with the 0(3) triplet as 
9 -> tt jij. 

Since the WZNW model is a conformal field theory, 
the 2-d 0(3) model at 9 w tt is a natural candidate for a 
slowly walking asymptotically free theory near a confor- 
mal fixed point [HI] . Following [3(3, EH , we define a run- 
ning coupling constant a(9,L) = g 2 (9,L) = m(9,L)L 
through the massgap m(9,L) in a periodic volume of 



spatial size L. At small L the coupling a(9, L) can be 
computed in perturbation theory. In this limit, it is in- 
dependent of the vacuum angle 9, which does not af- 
fect perturbation theory, and it agrees with the standard 
asymptotically free coupling constant of the 2-d 0(3) 
model. For large L, on the other hand, the coupling 
a.(0, L) is (9-dependent and can only be computed non- 
perturbatively. The corresponding /3-function is given by 
0(6, a) = —L8lo(9, L). Thanks to an ingenious use of 
the thermodynamic Bethe ansatz the finite- volume mass- 
gap is known analytically, both at 6 = [HI and at 
9 = 7r [33[. Although the 2-d 0(3) model is not inte- 
grable for other values of 9, by expanding around 9 = tt 
some interestin g an alytic results have been obtained even 
in that regime [29(. Due to the non- vanishing massgap 
M at 9 = 0, for large L the coupling o(0,L) -)■ ML 
increases linearly with the volume. At 9 = n, on the 
other hand, the massgap vanishes and the coupling ap- 
proaches a fixed point a(n, L) — > a* — n as one increases 
L — >• oo. In addition to the Gaussian fixed point at 
a = 0, the /3-function has another zero at a* = n, i.e. 
(3(9 = it, a = 7r) = 36]. Since the physics is sym- 
metric around 9 — tt, at the fixed point the /3-function 
just touches, but does not cross zero. The double zero of 
the /3-function thus corresponds to the lower edge of the 
conformal window. The parabolic form of the /3-function 
near this fixed point causes large logarithmic finite-size 
corrections: 
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These logarithmic corrections lead to a very slow ap- 
proach to the conformal fixed point. In our model, they 
are associated with a marginally irrelevant operator that 
breaks the 0(4) symmetry down to 0(3) [3J-|37[. 

Besides the running or walking coupling, the 9- 
dependence of the infinite- volume massgap is also of in- 
terest. Near the fixed point at 9 = tt it is given by [36( 

-1/2 



m(6, L -> oo) ~ |6»-7r| 2 / 3 |log(|6i-7r| 



(7) 



The determination of m(9, L) and a(9, L) requires non- 
perturbative calculations. These can be performed from 
first principles using lattice field theory. Simulating the 
2-d lattice 0(3) model at large 6 is very challenging due 
to a severe sign problem. Fortunately, the meron-cluster 
algorithm [38|, which is based on the Wolff cluster al- 
gorithm [39( , substantially reduces the sign problem and 
thus makes a numerical study feasible. For technical rea- 
sons related to the meron-cluster algorithm, we consider 
the 2-d 0(3) model on a triangular lattice of spatial ex- 
tent L and Euclidean time extent (3 > L as illustrated 
in Figure 1. The action is defined on nearest- neighbor 
bonds (xy), and is given by 



S[e\ = s(e x ,e y ), s(e x ,e y ) = \(l - 



■x ^-y 



), (8) 



3 



for e x ■ e y > — h and s(e x ,e y ) — oo otherwise. This 
action eliminates field configurations for which the angle 
between neighboring spins exceeds 120 degrees, which is 
essential for the success of the meron-cluster algorithm. 
The geometric topological charge density q( xyz ) <= hi|] 
[iol ] associated with a triangle (xyz) is given by 

Rexp{2-Kiq( X y Z )) = I + e x ■ e v + e v ■ e z + e z ■ e x 

+ ie x ■ (e y x e z ), R>0. (9) 

Here ^itq( xyz ) is the oriented area of the spherical triangle 
on S 2 defined by the three unit- vectors e x , e y , and e z . 
By construction, the geometric lattice topological charge 
Q[e\ = J2(xyz) 1(xyz) i s an integer. In order to determine 
the massgap, we consider the operator E(t) = J2 Xl 
where the sum extends over all points x = (x\,t) in a 
time-slice, and we define the 2-point function 



de x £7(ti)-£(i 2 ) 



x exp(-S[e] + i9Q[e\) ~ exp(-m(0, L)(t 2 - h)), 
Z{9) =]J f de x exp(-S[e] + t9Q[e\), (10) 

which decays exponentially with the 9- and L-dependent 
massgap m(6,L) at large Euclidean time separations. 

Like in the Wolff cluster algorithm, in the meron- 
cluster algorithm one first chooses a reflection plane per- 
pendicular to a randomly selected unit-vector r G S 2 . 
In a given update step, the spins e x are either left un- 
changed, or they are reflected to e x = e x — 2(e x -r)r. Spins 
belonging to a common cluster are reflected collectively 
[39| . In the meron-cluster algorithm, each cluster C con- 
tributes an integer or half-integer to the total topological 
charge Q. Thanks to the 120 degrees angle constraint, 
the cluster topological charge Qc is well-defined and does 
not depend on whether other clusters are flipped or not. 
Clusters with a half-integer topological charge are called 
meron-clusters, because they represent half-instantons. 
Flipping a meron-cluster changes the topological charge 
by an odd integer. At 9 = n where exp(i9Q) = (— l)^, 
this implies an exact cancellation between two contri- 
butions to the partition function. Taking this cancel- 
lation into account analytically rather than statistically 
via Monte Carlo leads to an improved estimator that ex- 
ponentially reduces the sign problem. In a multi-cluster 
algorithm, for general values of 9 an improved estimator 
for the partition function is given by 



z(0) 

Z(0) 



(Qcos(tfQ c )). 



(11) 



Similarly, for the 2-point function one obtains 

(e x ■ e y exp(i9Q)} = ((e x ■ r)(e y ■ r) ]J cos(9Q c ) 

c 

x [5^,c v + (l-6c m ,c v )tM(0Qc m )tan.{0Qc v )]}, (12) 
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FIG. 2. (top) The running coupling constant a[9,L) = 
m(9, L)L as a function of the scale set by the spatial size L for 
9/n = 0, 0.25, 0.50, 0.75, 0.92, 0.94, 0.96, 0.98 and 1. At 9 tt 
the coupling increases linearly as L — > 00, while at 9 — n 
it approaches the fixed point value tt (indicated by the hori- 
zontal line). At 9 w 7r the coupling is slowly walking near 
the conformal fixed point. The curves for < 9 < 3tt/4 show 
the effect of leading finite-size corrections exp(-mL) / %/ mL to 
the infinite-size massgap m. (bottom) The corresponding /3- 
function obtained by differentiation of our numerical results, 
supplemented by exact values at 9 = n from fH/. 



where the cluster C x contains x and C y contains y. 

We have used the meron-cluster algorithm to deter- 
mine the massgap m(9, L) for a large variety of 9- values 
and volumes L ranging from 6 to 100. By using at least 
four values of the bare coupling g, all results have been 
reliably extrapolated to the continuum limit. The cor- 
responding results for the running coupling a(9,L) are 
illustrated at the top of Figure 2. Within error bars, 
they agree with values obtained from the exact massgap 
both at 6 = and at 9 = tt [33]. While a(9,L) in- 
creases linearly for large L when 9 ^ tt, it flattens off for 
9 = tt. As anticipated, due to large logarithmic correc- 
tions the approach to the fixed point is very slow. For 
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9 » 7T, the coupling walks almost as slowly as at = 7r 
up to some distance scale, at which it starts running off 
into the linearly rising regime. While in a lattice con- 
text it is most natural to use a step scaling function 31| , 
here we prefer to discuss the /3-function that is familiar 
from the continuum, although its computation requires 
a spline-interpolation of the lattice data. The walking 
versus running of the coupling manifests itself in the /3- 
f unction shown at the bottom of Figure 2. For 9 « tt it 
walks towards the fixed point /3(9 — tt, a = tt) = before 
running off (linearly) to negative values. 

Figure 3 shows the finite-size scaling behavior 
of m(9, L)L + [tt — m(7r, L)L] as a function of 
MLt 2 / 3, 1 ' yj\ log(i/i )|, where t is the reduced coupling 
t = 1 — 6/tt, and to — 70. Large logarithmic corrections 
are removed in this difference between the finite- volume 
massgaps, so that the L —> oo value tt is enforced at the 
critical coupling t = 0. The fact that all data fall on 
a universal curve confirms the behavior of the massgap 
eq.([7]), and in particular the critical exponent 2/3 of the 
WZNW model. 

It is interesting to compare the behavior of the 2-d 0(3) 
model near the conformal fixed point at 9 = tt with the 
anticipated behavior of walking 4-d non-Abelian tech- 
nicolor gauge theories near the conformal window [25$ . 
Both theories are asymptotically free and conformality is 
thus limited to scales far below Ajjg- (i) First of all, in 
the 2-d 0(3) model the parameter that determines the 
distance to the conformal fixed point is the continuously 
varying vacuum angle 9, which does not get renormal- 



ized 24j, |25| , and which does not affect the /3-function in 
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FIG. 3. Finite-size scaling of m(9, L)L near 9 = tt as a 
function of MLt 2/3 /^/\ log(t/t )\, with t = 1 - 6»/tt, £ 
70 System sizes L have been chosen as multiples ofro/M, 
with ro = 0.26715356 JI2] . The data have been shifted by 
[tt — m(7T,L)L], to eliminate the large logarithmic corrections 
and enforce m(9 = n,L)L — tt. They fall on a universal 
curve, confirming the critical exponent 2/3 of the WZNW 
model in eq.^. This curve shows both the slow walking near 
n, and the linearly rising regime as L — > oo. 



the perturbative regime. In walking technicolor theories, 
on the other hand, the corresponding parameter is the 
discrete number of techniquark flavors or the size of the 
technifermion representation, which do affect the pertur- 
bative /3-function. Unlike in the 2-d 0(3) model, due to 
rcnormalization, one can then not directly compare phys- 
ical quantities between theories in and outside the confor- 
mal window, (ii) While large, logarithmic finite-size ef- 
fects are a characteristic of walking theories near the edge 
of the conformal window, as shown by eq.©, their ori- 
gin may differ. As we have seen, in the 2-d 0(3) model a 
marginally irrelevant operator breaks the enhanced 0(4) 
symmetry in the low-energy sector. In a 4-d non-Abelian 
technicolor gauge theory, it is not clear whether loga- 
rithmic corrections could come from a similar symmetry 
enhancement, or from conformal symmetry itself, (in) 
In the 2-d 0(3) model the 0(4) symmetry enhancement 
causes an 0(3) singlet to become light as 9 — > tt, in ad- 
dition to the 0(3) triplet, all having masses much below 
Ajt q, before both objects become massless "unparticles" 
41] at 9 = tt. It has been argued that particles with a 



mass much below Ajjg should also arise in walking tech- 
nicolor theories near the edge of the conformal window 
42 , 43j . These so-called technidilatons have been identi- 



fied with pseudo-Nambu-Goldstone bosons of a sponta- 
neously broken conformal invariance, which is still weakly 
explicitly broken by the scale anomaly at Ajfg. In the 2-d 
0(3) model, due to the Mermin- Wagner theorem, con- 
formal invariance cannot break spontaneously and thus 
the light 0(3) triplet and singlet are not expected to be 
pseudo-Nambu-Goldstone bosons. In particular, they are 
exactly massless at 9 = tt, despite the fact that confor- 
mal symmetry exists only in the low-energy sector and 
is explicitly broken at the scale Ajjg- (iv) In addition 
to the large finite-size effects, cut-off effects due to a fi- 
nite lattice spacing a may also be important. In the 2-d 
0(3) model, large logarithmic corrections to the expected 
0(a 2 ) effects mimic 0(a) behavior 44]. Similarly, in 
lattice investigations of technicolor gauge theories, one 
must control these cut-off effects. When one uses Wil- 
son fermions, which are theoretically cleaner than stag- 
gered fermions, without Symanzik improvement lattice 
artifacts are of order a. (v) Finally, in a non-Abelian 
gauge theory the Schrodinger functional [451 ] provides 
a definition of a running or walking coupling constant, 
which naturally replaces the coupling based on the finite- 
volume massgap that we use in the 2-d 0(3) model. 

In the end, our study of the 2-d 0(3) model near 9 = tt 
demonstrates that slow walking can indeed be studied 
accurately using Monte Carlo simulations, provided that 
lattice artifacts and finite- volume effects, which may both 
be large, are well understood and under good numeri- 
cal control. Besides technicolor gauge theories, it would 
be interesting to also investigate other models, e.g. 4-d 
Yang-Mills theories at non-zero 9, in order to further in- 
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vestigate the neighborhood of conformal fixed points. 
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